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Spontaneous and stimulated emission of an electron in the ion channel is studied. The emission 
processes are studied in the regime of high harmonic generation when the parameter of plasma 
wiggler strength is large. Like for conventional free electron laser, a synchrotron-like broadband 
spectrum is generated in this regime. The asymptotic expression for the radiation spectrum of the 
spontaneous emission is derived. The radiation spectrum emitted from axisymmetric monoenergetic 
electron beam is analyzed. The gain of ion-channel synchrotron-radiation laser is calculated. Use 
of laser-produced ion channel for efficient X-ray generation is discussed. 

I. INTRODUCTION 

An electron dynamics in plasma- focusing channel has important applications to new plasma technologies, such as 
advanced accelerators novel radiation sources, new types of lens jj. It is a key phenomenon for ion-channel laser 
(ICL) ion-ripple laser plasma- wiggler free electron laser(FEL) which are perspective candidates for high 
brightness X-ray radiation sources. Such radiation sources are strongly needed for experimental research in physics, 
chemistry, biology and in engineering . 

Resent experiments that explore the interaction of intense 28.5-GeV electron beam with plasma at Stanford Linear 
Accelerator Center (SLAC) 0, Q have shown that ion channel can be successfully used to produce broadband X-ray 
radiation. Moreover, the high density of the ions in the channel provides much higher wiggler strength than that 
provided by a conventional magnet wiggler. This leads to a more effective generation of X-ray radiation than in 
conventional light sources and could be used for the development of next generation of radiation sources. 

To create ion channel, electron beam have to interact with plasma in blow-out regime when the electron beam 
density, n;,, is higher than plasma density, Up. In this regime the electron beam charge quasistatically expels slow 
plasma electrons inside and around the beam volume and the ion channel is formed. Note that relativistic electrons of 
the beam are not expelled from the channel because of relativistic compensation of the beam electron charge force by 
the self-magnetic force |lOj| . The channel radius, ~ TQ^rib/np is much more than electron beam radius, tq, for dense 

nb 3> Up and narrow fc^ro ^ 1 electron beam where kp — c/ujp is the plasma skin depth, ujp — {^inripe'^ /m)^^"^ 
is the plasma frequency, e is the electron charge, m is the electron mass and c is the speed of light. If all plasma 
electrons are expelled from the channel then the restoring force on the beam electrons due to the ion charge is given 
by Gauss's low and in cylindrical geometry is: 

Fres = (1) 

where rj^ is the vector from an electron to the channel axis. The beam electrons in the ion channel will undergo betatron 
oscillations caused by this force. The wavelength for small betatron oscillations is close to Xb = 2-K/kb ~ ^-K^flFf jkp, 
where 7 is the relativistic factor of the electron beam 12]. 

It is well known that accelerated charges emit electromagnetic (EM) radiation [l^l ■ Therefore the electrons under- 
going betatron oscillations in ion channel will emit short-wavelength EM radiation. Some features of this radiation 
spectrum have been studied in Ref. 0jllil3- The wavelength of the radiation is close to A ~ Ab/ (27^) for small- 
amplitude near-axis betatron oscillations. If the amplitude of the betatron oscillations becomes large then electron 
radiates high harmonics. If the plasma wiggler strength, 



K = -jkbTQ = 1.33 • 10~^°\/7"e [cTO-3]ro \\xm\ , (2) 

is so high that K ^ 1 then radiation spectrum becomes quasi-continuous broadband, where tq is amplitude of electron 
betatron oscillation in the channel. The frequency dependence of the radiation spectrum becomes similar to the one 
of synchrotron radiation spectrum which is determined by the universal function 



S (cj/ujc) = i^^/uic) K^/3{x)dx, (3) 
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FIG. 1: The synchrotron radiation function S{x) versus x. 



where ujc is the critical frequency (see Fig. For frequencies well below the critical frequency {uj <C uJc) the 

spectrum increases with frequency as w^^'^, reaches a maximum at ^ 0.29wc, and then drops exponentially to zero 
above ljc- The critical frequency for a relativistic electron in an ion channel is |14| 
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Wc = -^l^cr^kl = 5.2 • 10~^**7^ne [cm"^] tq [^Jim] keV. (4) 

Because of the strongly relativistic motion of the electron the emitted radiation is confined in very narrow angle 
6 ~ K/^. Synchrotron radiation in ion channel has been observed in recent experiments^ 
The averaged total power radiated by an electron undergoing betatron oscillations is 

2 

(Ptotai) ^ \^m^Xrl (5) 

where A'';, is the number of the betatron oscillations performed by the electron. We can introduce also the stopping 
power of an electron, that is the energy loss of an electron per unit distance 

9 MpV 

Q = {Ptotai) /c ^ 1.5 • 10-45 (^^^ [^^-3] _e^^_ 

cm 

The averaged number of the photons with averaged energy hujc is 

(Nph) ^ ^^NbK ~ 1.02 • 10-^ NbK. (7) 
9 he 

It follows from Eq. Q that the radiatedpower is proportional to the squared density of ions in the channel. This fact 
has been confirmed in the experiments T|. As it was mentioned above the ion density in the channel have to be less 
than electron density in the beam. This leads to the serious limits on the gain in radiated power. One of the ways 
to overcome the limits is to use laser-produced ion channel. The ion density in such channel can be about lO^^cm"'^ 
that is in 5 order higher than ion density in the channel produced in beam - plasma interaction experiments 
So use of ion channel produced by laser pulse could increase the power of X-ray radiation in 10^" times! More 
detailed discussion on use of relativistic laser channel for X-ray generation will be p resented in Conclusion. 

Spontaneous emission in an ion channel has been studied in detail in Ref. |l4| . The general expression for the 
spectrum has been derived. It is a complex expression that involves the sum of products of the Bessel function. 
Numerical evaluation of the spectrum becomes difficult in the limit if 3> 1. The simple asymptotic expression for the 
angular dependence of the radiated spectrum has been derived for this limit only for directions that are perpendicular 
to the plane of the betatron oscillation. The one of the purposes of our paper is to calculated the total angular 
dependence of the radiated spectrum. 

Resonance interaction between EM radiation and the betatron oscillations of electron beam in ion channel leads to 
the bunching of the electron beam and then to the amplification (or damping) of EM radiation. It is a stimulated 
emission (or absorption). Stimulated emission is a basic process in ICLs 3] and FELs [17|. Reverse process - 
stimulated absorption leads to the direct laser acceleration 18] and to the magnetic field generation in relativistic 
laser channel. Unfortunately stimulated emission (absorption) in ion channel has not yet been explored in the limit 

^ 1 that is another objective of our paper. 
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The paper is organized as follows. In Section II the motion of an electron is studied and Hamiltonian formulation of 
the problem is presented. In Section III, the spontaneous emission from electrons undergoing betatron oscillations in 
ion channel is analyzed. The general asymptotic expression of the radiation spectrum for arbitrary angular dependence 
is derived for K ^ 1. The spectrum averaged over azimuthal angle is calculated for axially symmetrical electron beam 
in the limit K ^ 1. Section IV discusses the stimulated emission processes. The gain of ion - channel synchrotron - 
radiation laser is derived. A summary discussion is presented in Sec. V. 



II. ELECTRON DYNAMICS IN ION CHANNEL 



Relativistic equation of electron motion in cylindrical ion channel is 

^=F_, (8) 

where F^es is the restoring force defined by Eq. It follows from Eq. © that momentum along channel axis Pz 
is a constant of motion. First we will consider the radial betatron oscillation as it takes place when the center of an 
electron bunch moves along the channel axis. Assuming that py — we get the equation for x -coordinate: 

where we introduce the constant of motion 7^ ~ ^+pt E^nd use the dimensionless units, normalizing the time to UJ~^, 
the length to c/ujp, the momentum to mc. 

As Hamiltonian does not depend on time it is another constant of motion 

2 2 
H = j+^= const = 7. + ^, (10) 

where tq is the amplitude of the betatron oscillation. We can express 7 as function of x from the obtained relation to 
resolve Eq. . Then transversal motion of electron can be reduced to the oscillations in effective potential 

U (x) = ^ 2 . (11) 

The oscillations can be described in implicit form as follows 



Vj^^ + 4 ^ / .fx 
t — tq E arc sm — 



^ ^^(arcsin(-),^:— ), (12) 



i^y/iy'^ + 4 V V ''O 

where E{x, k) and F{x, k) are the Elliptic integrals of the first and the second kinds 2(3|, respectively, and u is 

- r2/27,. (13) 

The period of betatron oscillations is 



2^ 

Tfc = — = 4ro 



(14) 



where K{x) and E{x) are the complete Elliptic integrals of the first and the second kinds |23|, respectively. In the 
limit ^2^1, parameter v is the ratio between the longitudinal and transversal energy of the electron 

^Plhl^p\lpl, (15) 

where p^ is the maximum of the transversal momentum which is at the channel axis a; = 0. In most applications 
transversal moment of the electron is much less than longitudinal one, so we assume that <C 1. Then we can use 
expansion in v to describe betatron oscillations: 

(16) 



/27I 
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X (t) — tq sin (ujbt) — ^'OTTTi'^ sin (3ujbt) + 
64 



(17) 



Using E g. 11 OH we obtain the relations for the electron orbit in the zeroth order in v that coincides with ones calculated 
in Ref. ^ 



1 



(18) 



X (t) ~ To sin {ujbt) , 



(19) 



y {t) ^ 0, 



(20) 



z (t) ~ Zo + 



1 - — \t- To — - sm {2ujbt) . 
4 / 7z 8 



(21) 



Notice that parameter coincides with the expression kbro introduced in Ref. [T^ and plasma wiggler strength 
parameter can be expressed through v as K — p±. 

More general regime of the betatron motion when Py and an electron orbit is not plane has been considered 
and classified in Ref. [l9l | in the limit p^ 3> p± ■ Equations of motion take a form in this case 



x{t) 



Px 



sin (wfji) , 



(22) 



vit) 



Py 



sin {uJbt + ip) 



(23) 



Z (t) ~ Zo + 



Pi 



p 

87, 



2 sin {2ujbt) 



87.2 



^ sin {2<jjbt + 2tP) 



(24) 



where '0 is the phase difference between oscillations along a;-axis and y-axis, p^ and Py are the maximum of the 
electron momentum along cc-axis and along y-axis, respectively, that occurs at the channel axis (cc = 0, y = 0). If 
angular momentum of the electron, L — pyX —PxV is equal to zero, then electron executes radial harmonic oscillations 
through the origin with amplitude tq = 2^/H. If L ^ ±imax = ^H/cub then an electron performs circular motion 
with radius vq. In the general case (an arbitrary value of — imax < L < Lmax) the electron trajectory is an ellipse-like 



and is confined between the maximal radius, ^ IH + y/H^^^ujfL^ ] , and minimal radius, w I iJ — \JlP 



III. SPONTANEOUS EMISSION 

Using Eqs. H22I) - H24I) for electron trajectory, the energy spectrum radiated by an electron can be calculated p^[T3|. 
The total radiation flux can be separated in two independent components with polarization in the ee and directions, 
where unit vectors eg and correspond to spherical coordinate system: x — r sin cost/), y — rsin^sin^, z = rcosO. 
Then energy radiated per frequency oj per solid angle in the direction k = (w/c) (e^, sin 9 cos (p + ey sin ^ sin -|- cos 9) 
during the interaction time T is jl4| : 

dWspon,j _ ^'.^2|j^.|2^ (25) 



T/2 



' dx dy dz 

Ig = / I — COS COS (/) + — COS 6* sin — sin^? ) e'*c?i, (26) 

J \dt dt dt ' 

-T/2 
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FIG. 2: Schematic of synchrotron radiation in an ion channel. Open circles show the points on the electron trajectory when 
the electron emits in the direction of k. 

T/2 

I^^ I (^^sin0-^cos(/))e'*di, (27) 

-T/2 

* = 7/ [i - x (t) sin 9 cos 4>-y (t) sin 6* sin - z {t) cos 9], (28) 

where j = 9,<j) is the polarization index, the electron trajectory is given by Eqs. H22|l - 124|l and 77 = w/wp. The final 
results can be expressed as double infinite series of the Bessel function products (see Eqs. (32) - (41) in Ref. \lM or 
can be expressed by the infinite series of the generalized Bessel function introduced in quantum electrodynamics 21]). 
Unfortunately the series converge slowly in the limit if ^ 1 that makes the numerical evaluation of the radiation 
spectrum difficult. Energy spectrum and angular dependence of radiation have been derived only for directions that 
are perpendicular to the plane of the betatron oscillation (i. e., for (f) ~ 7r/2) To extend this result we will use 
saddle point method to evaluate integrals (^5)1 and (|?7l) . 

It is well known that radiation of accelerated relativistic electron is beamed in a very narrow cone in the direction 
of the electron momentum vector, and is seen by the observer as a short pulse of radiation as the searchlight beam 
sweeps across the observation point [T^ (see Fig. So the time moments when the electron momentum, p, is 
directed along wave number, k, give the main contribution in integrals ( and (j?7|) . As the pulse duration is very 
short it is necessary to know the electron momentum and the electron position over only a small arc of the trajectory 
whose tangent points in the direction that is close to the direction of k. Then we can expand integrand in Eqs. I|2t)|) 
and H27|l about this time moments and perform integration. This approach implies that we approximate the part 
of the electron trajectory near these time moments by the arc of a circular path |l3l| . In this case the radiation 
will have well-known synchrotron-like spectrum. It is noted in Ref. [23| that another necessary (the first is that the 
electron have to be relativistic) condition for use of the synchrotron radiation approach is that the electron deflection 
angle should be much more than simultaneous angle spread to which radiation is emitted. The electron momentum 
oscillates in cone angle ^ V^IVz (electron deflection angle). The radiation of the relativistic particle is confined to 
angle I/7 |l3l |. So the validity condition for synchrotron radiation approach is Pi_/pz 3> I/7 that is equivalent to the 
condition for high harmonics generation discussed in Introduction K ^ 1. Therefore the radiation spectrum can be 
approximated by the synchrotron one if 

7>P±>1. (29) 

This condition is easily satisfied in experiments. 

The arguments presented above justify the use of the saddle point method [23| to evaluate integrals (|26|l and l|27|l . 
We can expand phase ^E" about the moment of time f„ = LUbtn to the third order: 

* = *o + bi{(, - e.) + h2{£, - £.n? + - in)\ (30) 

*o = * (^„) = aoC« - axsin(^„) - sin(^„ -I- V") 

+azx sin (2^„) + a^y sin (2^„ -I- 2?/^) , (31) 
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h = 



= ao- ax cos (^„) - ay cos + tp) 



+2aza: COS (2^„) + 20;^^ COS (2^„ + 2-0) , 



(32) 



ax sin (^„) + sin (^„ + ■0) 



-4az2; sin (2^„) - 40^^ sin (2^„ + 2-0) , 



(33) 



COS (^„) + ay cos (^„ + -0) 



-Sa^a; cos (2^„) - 8a^y cos (2^„ + 20) 



(34) 



where 



ao 



UJb 



Pz r pI+pI . . 

— 1 . „ sm u cos ( 

7. \ 



(35) 



V Px . a 

a^ = sm a cos < 

W6 7^ 



(36) 



a„ = sm y cos < 

Wfc 7^ 



(37) 



^ Px 
Jz 872 



COS 9, 



(38) 



V Pz Py a 
^zy = —^COS0. 



(39) 



Then we can expand the pre-exponent factors in Eqs. H26|) and (|27|l about the moment of time ^„ = LOhtn o the first 
order: 



dx 



dy 



dz 



COS0COS0 + — cos sin — — sin0 = Bg^n + Dg^niS, — Cn) 



dt 



dt 



(40) 



dx 

in 



dy 

I 

dt 



i0 = B^,n + D^^ri{i - in), 



(41) 



where 



— cos in COS H — - COS (^„ + tp) s'm ) COS 6 

Iz Iz 



— smf 



.Pz 



pI 



PZ 



2 cos (2^„) - ^ cos (2e„ + 20) 



(42) 



— sin in cos + — sin (^n + 0') sin 

7z 7^ 



COS( 



Pz pI 



Pz PI 



£1 ^ sin sin {2in)-—ir^ sin sin (2C„ + 20.) 
7z 272 7^ 27^ 



5<A,n = — sin cos (in)-— COS COS (^„ + 0j) . 

7^ 7^ 



(43) 
(44) 



7 



sin sin + — cos (j) sin 



7^ 



7z 



(45) 



Notice that it is sufRcient to keep the leading term in the pre-exponent factor while only terms that are much less 
than unity can be neglected in the exponent argument. 

The main contribution to the integral comes from the neighborhood of the saddle points specified by d^'/d^ — 
[2^. The first-order term in phase expansion (|30|l can be written as follows 



UJb 



1 



kp 



2ujb \r 



1 



(46) 



where tp is the angle between k and p. It follows from Eq. 1)46(1 that d'i'/d^ is minimal and close to zero a.t (p — 
when the electron momentum is directed along k that agrees with the qualitative argument presented above. 

For simplicity we assume that Py = that is the electron orbit is plane and the betatron oscillations is radial. It 
follows from Eq. that the values of ^ whose neighborhood gives the main contribution to the integral are defined 
by the relation 



cos^„ 



■ tan 9 cos i 



Px 



(47) 



It is seen from Eq. (|47|l and Fig. 3 that the number of saddle points is 2Ni, — uobT/n that is the number of times when 
the direction of the electron momentum and the direction of the wave number coincides. It can be shown that 
the second-order term can be neglected in Eq. ((30(1 . Then Eas. ((26(l and 1(27(1 take a form: 



Ij = — exp (z^o,n) Rj.n, 

71—1 

/ + 00 
dSn {Bj^n + Dj.nSn) exp [i6l,„S„ -I- 163, nS^] 
-00 



(48) 
(49) 



where j = 9,(f> is the polarization index, Sn = £,n — Co, ^o,n is the value of phase in the n-th saddle point. 

It follows from the definitions of saddle point Eq. ((47(1 and from Eqs. II3U(I - l(45() that Rj n ~ Rj,n-i = Rj, 
Bj n = Bj_n-i = Bj, Dj n = Dj n^i = Dj, bm,7i — &m,ji-i = • Performing integration in Eq. 1(49(1 we obtain 



(50) 




where Ki/3{x) and -^^2/3(2;) are the modified Bessel function |20j. In the synchrotron regime of radiation lu 
and ^o,n ^ 1 then we can write for large number of the betatron periods {Ni, :s> 1 ) 



p 



2Ni 

^ exp (j^'o.n) 
n=l 

Using Eqs. (^HJ, lEOI) and condition (EHJl we finally get 



2Nh. 



spon.O 

diodfl 



2N, 



9 9 2 



[sin 9 siv? (t>Ki/3 (q) + cos 0if 2/3(9)] ^ , 



(51) 



(52) 



dWspon,4, ^^'^V^Xsin^f/) |- . at^ f \ /-r^ f \^'^ /t,o\ 

^ -2Nb — [sm9 cos (t)K^i:i{q) - ^K2/3[q)\ , (53) 



where 



(54) 



V Px/lz ~ sin 9 cos^ (/) 

is the curvature radius of the circular path that is used to approximate the part of the electron trajectory where 
electron emits in the direction of k, 

X = 7^^ + sin^ 9 sin^ 0, (55) 
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FIG. 3: Angular distribution of tlie spontaneous emission spectrum ^J^°" (arbitrary units) versus angles Ox and dy from a 
single electron with 7^ — 1000, K = px = 20 for u — O.SuJc- 



3/2 



The total radiation of the spontaneous emission from the electron in the channel is 



(56) 



spon 



dW. 



spon, 6 



dW. 



spon,(p 



dujdVL 



duodVL 



dujdQ 



2N, 



sin^ 9 sin^ 1 



X 



(57) 



This is the general expression for the angular distribution of the radiation emitted by an a relativistic electron in ion 
channel and this is the one of the main results of the paper. 

Let us consider some limiting cases. It follows from condition H29|l that 0^1. Then in the limit (f) ^ n/2 Eq. H57() 
takes a form 



dWspon _ ^ 6e2 
duidVt 



9 9 



(1+7202) 



7, 



,2/12 



■^1/3(9) +^2/3(9) 



(58) 



that coincides with the asymptotic spectrum emitted by the relativistic electron in the channel for (p — ■n/2 (see 
Eq. (64) in Ref. 0|). It was discussed above that in the limit K ^ 1 the radiation emitted from an electron at the 
given moment of time is similar to the synchrotron radiation emitted from an electron in instantaneously circular 
motion with the same curvature radius. Really, introducing notation Lp — sin sin and using relation 7 ~ 7^ ~ 
we can reduce Eq. H57(l to the known form 



1 

7 



^2/3(9) + 



1/72 + 



(59) 



that coincides with the expression for energy radiated by a relativistic electron in instantaneously circular motion 
with radius p per unit frequency per unit solid angle 27TCOS(pd(p after 2Ni, revolutions (see Eq. (14.83) in Ref. 13]). 

To visualize our results we use new angle coordinates 9^ = sin cos 0, 9y = sin sin instead of the spherical one 
{9, If) (see Fig. Eland Fig. OJ. It is seen from Fig. |3| that there is no radiation for 0^ > 6'max because the argument 
of the Bessel function in Eq. (|57() goes to infinity and the Bessel function goes to zero for 0^ = 0max- Therefore 
in our approximation the emission angle along x-axis is confined to the electron deflection angle 0max — Px/Pz- 
However it is evident that electrons with maximal deflection angle radiation Px/pz emit radiation up to the angle 
6* = ± (6'max + 7~^)- The emission angle is confined to the angle ^ I/7 in the direction of y-axis which is normal to 
the electron orbit plane. Hence our results agree with the qualitative analysis in Ref. p^. 

Averaging Eq. H57|) over 9y we obtain 



-^^d9y - —N, 




(60) 
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where S (x) is the universal function mentioned in Introduction. Expression ()60|) can be considered as a radiation 
power from electron flow in which electrons are evenly distributed over y-axis (over direction which is normal to the 
betatron oscillation plane). 

Let us now consider a monoenergetic, axisymmetric electron beam. Radiation spectrum from electron beam with 
electron distribution function / (p) is defined by the relation: 

^ *yspon 



diodQ 



Let all electrons have the same longitudinal and transversal energy before interaction and the electron distribution 
function is 



f{p) = S{p,^^,)6{^pl+pl-p^). (62) 
Then the radiation spectrum from the beam takes a form 



/ dWspon \ _ f^"" d(j) dWspon ^g^^ 



In the limit of near-axis radiation 9 ^ l/7z the radiation spectrum is 



fe) (fr^) -V. (f?^) ^ 

In the reverse limit 6 ^ l/7z the main contribution to Eq. (|63|l is given by the small values of <C 1 that correspond 
to the minimum of q. The radiation spectrum in this limit is 

/ dW^pon \^,V3 e^7. V277 \ 

\ dLodn / ~ 2 TTC9 ^3 ^3/2^ ^2 _^2g2 J • ^ 

Note that the obtained expression coincides with Eq. (|^ . 



IV. STIMULATED SYNCHROTRON RADIATION IN ION CHANNEL 



When EM wave of approximately the same frequency as the spontaneous emission is driven through the ion channel 
simultaneously with the electron beam, a significant exchange of energy between the beam and the wave can occur 
and can lead to the coherent efficient amplification of the wave energy. This amplification can be explained in terms 
of the stimulated emission processes which is determined the operation of ICL. Unlike ICL theory Q, here we will 
consider regime of strong wiggler {K 3> 1) when the emission process is close to synchrotron one. This is the regime 
of ion-channel synchrotron-radiation laser (ICSRL). 

The difference between spontaneous and stimulated emissions is the following. The radiation fields generated by 
an electron undergoing betatron oscillations has a phase which depends on time of arrival of the electrons at the 
channel entrance. The fields produced by different electrons in a uniform input beam have a random phase relation 
to each other and sum up incoherently. This lead to incoherent spontaneous radiation. Contrary to the spontaneous 
emission, electrons in ion channel can be driven by an external EM wave in synchronous oscillation, the phase of which 
is no longer random but locked to the phase of the wave. As a result the external wave causes the bunching of the 
electron beam and more efficient interaction between the beam and the wave. Then the radiation fields of different 
electrons sum up coherently to each other and to the external wave, leading either to a decrease or to an increase of 
the power, depending on whether destructive or constructive interference is realized. So the external wave is either 
absorbed or amplified. In quantum approach the amplification/absorption can be viewed as a transition, forced by 
the external wave, between the quantum states of the electron in the ion channel with photon absorption/emission. 
An amplification process of this kind is called a stimulated emission. It is known from laser physics that stimulated 
emission can be much more efficient and powerful than spontaneous one. 

As known in quantum physics there is the relation between the spontaneous and stimulated emission. Using 
this fact the elementary quantum methods based on the Einstein's coefficients have been used to study the instability 
of EM waves in cosmic plasmas J23]. Particularly, the synchrotron instability in a cold magnetoactive plasmas has 
been identified 26] . The relation between spontaneous and stimulated emission of the electrons in undulators is called 
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Madey's theorem in the theory of FELs |23 ■ It can be considered as extension of Einstein's coefficient method to the 
classical limit. Generalized Madey's theorem ^3 enables us to reduce the problem of ICSRL gain to the problem 
solved in Sec. Ill that is the calculation of the power of spontaneous emission in ion channel. 

To use Madey's theorem we should formulate the problem within Hamiltonian approach. An electron motion in an 
ion channel and in EM wave can be described by a relativistic Hamiltonian 



H=^l + {p^+ A,y+pl + (66) 

where is the vector potential of the wave with 9 or polarizations (z — 6,(j)): 

Ag = Aq {ex cos9 cos(j> + ey cos9 sm(j> ~ ezSinO) 

X exp [irjt — ir] (k ■ r) /k] , (67) 



A^ = ^0 {^x sin 4> — ey cos (j)) exp [irjt — irj (k ■ r) /k] . (68) 

Hamiltonian H66|l is again written in the dimcnsionless units, normalizing the time to ujp^, the length to c/ojp, the 
momentum to mc, the vector potential to mc^ /e. As usual we assume that the time of arrival of the electrons at the 
channel entrance is random. Assuming that external EM wave is weak we can consider it as perturbation and use 
the perturbation theory to calculate the work done upon the electron beam by EM wave with j-polarization. Then it 
follows from the generalized Madey's theorem (see Eq. (5.52) in Ref. that this work per beam electron is 

where jij is the first-order work done upon a single electron moving along the unperturbed electron trajectory (t) 
by EM wave with j-polarization 

71.- j ^ dt^.A4i,r°(0], (70) 

7o and are the electron energy and the transversal momentum of the electron before the interaction. The unper- 
turbed electron trajectory r°(t) is determined by Eqs. H22|) - (|24|l . Averaging in Eq. H69|l implies the averaging over 
the time of arrival of the electrons at the channel entrance. Mathematical statement of the Madey's theorem is that 
the second-order quantity, {Wj), is proportional to the average squared first-order quantities, 71 j. Therefore Madey's 
theorem essentially simplifies calculations in the framework of the perturbation theory. Using instead of variables 70 
and px the initial value of the electron momentum p = (p^,pz), Eq. (|69|) can be also rewritten in more symmetric 
form 



Ik g(7L -) 

2 k' ap 



(71) 



Similar to the FEL theory we introduce the incremental gain of ICSRL as a ratio between a power generated by the 
electron beam and the incoming EM wave power 

where rib is the density of the electron beam, re = mf? je^ is the classical electron radius, A = 2itc/uj is wavelength of 
EM wave. It should be noted that to calculate {Wj) we consider the given EM wave and do not consider the dynamics 
of EM wave during the interaction. Hence our calculations is valid for F <C 1 . This regime of interaction is called in 
the FEL theory as a small-signal small-gain regime . 

It follows from Eqs. H26|) . H27|l and H70I) that 7i,j is proportional to Ij. That is the particular case of the general 
reciprocity relation between the far field of a moving electron and the work done by a plane EM wave on it |l7j . 
Therefore we can express quantity {li j) through the energy of the spontaneous emission energy radiated per frequency 
per solid angle 
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and can express the gain in term of 



^~ TT \4tt'^coj^J \k' dp) dujdn ' ^ ' 

It has been noted in previous Section that efficient spontaneous emission in direction k takes place only at a 
short moment of time when the electron moment is directed along k. Therefore we can conclude from Eq. H74() 
that interaction with EM wave propagating in direction of k is only possible at the same moments of time. It was 
mentioned in Sec. Ill (see Fig. 3) that the number of the interactions moments is 2^*";,. The number of the electron 
oscillations in EM wave between interaction moments is extremely large: Nemw — <^l^b — '^d^b — 37^i^/2. For 
example for SLAC experiments Nemw — 9.2 • 10^. So we can consider the electron phases as random before each 
interaction moments and, therefore, can consider each interaction moments independently. Using this fact Eq. 174|) 
can be rewritten as follows 

IT \ait^cui) \k dp)\ diodn nJ- ^^^^ 

For simplicity we assume that Py = 0. To take derivatives in Eq. H75|l we should present ''^^^"^'^ as a function of 
momentum, p, and use the simplifying assumptions Py = only after performing differentiation in Eq. H75|). This 
is because the electron motion in the ion channel and in EM wave is no longer plane. Although the unperturbed 
betatron oscillations is in the plane y = the action of EM wave leads to the small oscillations along y-axis in the 
first-order of perturbation theory. 

It follows from Eqs. (j2Sl, (EHl and (jnH) that 



f Uj'^\ ^ dWspon.j 1 _ 2 2|r>|2 



-2rm- (76) 



b 



Then using Eqs. H32|l - H50() we can perform differentiation in Eq. (|74l) . To do it we have to take into account that 7^, cjf, 
are the function oi Pzi 7z {pz) = ^/l+pl, cub (Pz) = ^/ y^lziPz)- After performing differentiation we can put py ^ 0. 
Then simplifying the obtained expression with help of MATHEMATICA the gain of ICSRL can be derived 

r, = ^^ V2AVn,/sin^^ [sin0 cos^K,/,iq) - 7^X2/3(9)] (77) 
X [VX {SqTo + sin e^Ti) K,/:,{q) 
- (Sq Sine cos 0To + x'/'Ta) ^^2/3(9) 

To = sin^ 6 cos^ </> [Spj^ — 8p^7z sin cos (f> — PxJz sin^ cos 

+27^ sin^ 6* cos (ji (1 + 2 cos^ 0)] - 3pl , 
Ti = 6pl cos (j) + 2pl-fz sin 9 {3- 8 cos^fj)) 

-2pj;7^ sin^ 9 cos^ (j> - 27^ sin^ 9 cos^ (/i(l - 4 cos^ 
T2 — 3p'^ — Spl^z sin 9 cos — Pxlz sin^ 9 cos^ (/> 

+27^ sin^ 6* cos 0( 1 + 2 cos^ </>) , 



2 



2^ „ „4 



Te = Nb ^^^J-""- [sin 9 sin^ (t>K^/M + Vx cos (pKyM] (78) 
Stt XPx 

X [Vx(3gcos(/.To+sin0V3eT3)ifi/3((?) 
+ (Sqsin^sin^ (t>To + X^'^T^ K2/z{q)\ , 

T3 = 6pl - IGpIjz sin 6* cos - 2pxj^ sin^ 9 cos^ 

+47^ sin^ 6* cos (1 + 2 cos^ (f) , 
T4 = 3pl cos (f> + 2pIjz sin 6* (3 - 4 cos^ 0) + p^Jz sin^ cos^ 

-47^ sin^ 9 cos^ cj) sin^ ^i. 
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FIG. 4: 

a) Te (arbitrary units), for ^-polarized EM wave with uj = luc versus angles 6x and 9y for electrons with 7^ — 500, K — px — 10. 

b) The domains of the angles 6^ and 9y where ^-polarized EM wave with lo = loc is amplified by the electrons Fg < 0, grey 
region) and where the EM wave is absorbed by the electrons (Fg > 0, white region) for 7^ = 500, K = = 10. The angles are 
given in radians. 

c) F^ (arbitrary units), for 0-polarized EM wave with co = Uc versus angles 9x and 6y for electrons with 7^ = 500, K = px = 10. 

d) The domains of the angles Ox and 9y where </)-polarized EM wave with uj = Uc is amplified by the electrons (F^ < 0, grey 
region) and where the EM wave is absorbed by the electrons (F^ > 0, white region) for 7^ = 500, K — px ~ 10. The angles 
are given in radians. 

We have also checked the obtained results by the numerical difFerentiation for some values of parameters. EM wave 
can be amplified only if the wave propagates in some small angle to the axis z. It is seen from Fig. ^ that there is no 
amplification of EM wave when the wave propagates along the channel axis. 

Eqs. H53|) and (|52|l have been derived under assumption that all electrons have the same momentum before inter- 
action. Now we will consider a monoenergetic, axisymmetric electron beam with electron distribution function given 
by Eq. (|^ . ICSRL gain for such electron beam is defined by the relation 

(r,)-^''^^r,(p,0). (79) 

We have performed integration in Eq. H79|l numerically. It is seen from Fig. [S] that there is no amplification of EM 
wave by the beam (j) for given parameters. We are not able to find the wave amplification by the axisymmetric electron 
beam at least for considered beam parameters. 

V. DISCUSSION AND CONCLUSIONS 

First we would like to discuss the use a laser-produced ion channel for X-ray generation. As it was mentioned in 
Introduction an ion channel in plasma can be produced by an electron beam itself. However in this case the plasma 
density have to be less than beam density. Unfortunately the density of relativistic beam cannot be very high because 
of technology reasons. This leads to the serious limitation on the gain in radiated power as the power is quadratic 
in plasma density. Use of high-power laser could overcome this limitation. High-power laser pulse can expel plasma 
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FIG. 5: 

a) (r^) (arbitrary units) for axisymmetric electron beam with 7^ = 500, p± = 10 and (/)-polarized EM wave versus angle and 
normolized frequency ui/uc. 

b) (Tg) (arbitrary units) for axisymmetric electron beam with 7^ = 500, p± = 10 and ^-polarized EM wave versus angle 8 and 
normolized frequency ui/uc- 




FIG. 6: Schematic of the spontaneous emission from external electron beam undergoing betatron oscillations in a laser-produced 
ion channel. 



electrons by ponderomotive force and create the ion channel behind the pulse Moreover, in strongly nonlinear 
"bubble" regime electrons are completely evacuated from the first half-plasma wave excited behind the laser 
pulse. The ion density in this "bubble" is higher in many order of magnitude than that in the ion channel formed 
in the beam - plasma interaction. For example, the ion density in the "bubble" can be as high as lO^^cm"^ [Tsl 
that is in 10^ times higher than that in the beam - plasma interaction experiments . Therefore the radiated power 
in laser-produced channel can be in 10^° times higher than that in the experiments for the same number of the 
betatron oscillations. The "bubble" moves with group velocity of the laser pulse which is close to the speed of light. 
Relativistic electron bunch injected into the "bubble" can propagate inside the "bubble" over a long distance. Hence, 
in spite of the small length of the "bubble" the electrons can oscillate in the "bubble" a lot of time. 

It should be noted that the self-generated forces of the electron bunch as well as of the relativistic plasma electrons 
in the column can be neglected because they cancel each other IQJ . The condition to neglect the self-generated forces 
is TT-e ^ n.ij'^ that is easily satisfied in modern experiments. The number of the betatron oscillations performed by an 
electron can be calculated as a ratio of the time during which electron passes through the "bubble" to the period of 
the betatron oscillation: 

2TTy/2-f{c ~ Vgr) 

where Vgr is the group velocity of the laser pulse. 

Some estimations for parameters in possible and carried out experiments are presented in Table I. To calculate 
parameters in Table 1 we assume that L ~ 30/im according to the numerical simulations jl^ . The group velocity 
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of the laser pulse in rarefied plasma (w^ ^ cUp) can be estimated as Vgr/c ~ 1 — w^/ (2w£), where is the laser 
frequency [ig. Then Eq. ^ can be rewritten as follows 

Al [M"^J V [cm •^J V Al / 

where Al is the laser wavelength. It is seen from Table I that use of the laser-produced ion channel can dramatically 



TABLE I: Parameters of SLAC experiments (a); parameters of interaction between 30-GeV electron beam and laser- 
produced ion channel (b); parameters of interaction between 30-MeV electron beam generated in laser-plasma interaction [l^ 
and laser-produced ion channel (c). 





ro 


7 


ni[cm ^] 


K 


LUclMeV] 


Nt 


^ L cm J 


Nph 


a 


40 


6- 10"* 


2 ■ 10" 


17 


0.12 


1.5 


2.5 ■ 10'** 


0.24 


b 


10 


6- lO"* 


1 ■ 10" 


1030 


1889 


2 


5.3 ■ lO"* 


21 


c 


10 


6- 10^ 


1 ■ 10" 


103 


0.76 


20 


5.3 


21 



increase the power of X-ray emission. 

In this paper we have studied spontaneous and stimulated emission from electrons undergoing betatron motion in 
ion channel. We calculate the period of nonlinear betatron oscillations. The method based on the Bessel function 
expansion is used in Ref. to calculate the spectrum of the spontaneous emission in ion channel. However, in 
synchrotron regime of emission [K >> 1) the angular distribution of the radiation has been derived in this paper only 
in the direction which is perpendicular to the electron orbit plane. We extended this result to the arbitrary directions. 
The generalized Madey's theorem was used to calculate the electron energy gain of ICSRL. The calculation shows 
that there is no wave amplification when the EM wave propagates along the channel axis and the amplification takes 
place only when the wave propagates at some small angle to the channel axis. 

Our numerical simulations show that there is no amplification for axisymmetric electron beam at least for considered 
parameters of the beam. However further analysis is required to justify this conclusion. It may be possible to overcome 
this difficulty by appropriately tailoring the electron beam, that is, a narrow electron beam could be injected off-axis 
such that all of the beam electrons execute approximately the same betatron orbit. 

To calculate the radiation spectrum and the gain of ICSRL for electron beam, we used a very simple distribution 
function. Further investigations should include the more realistic electron distribution functions. The gain of ICSRL 
was calculated in the small-signal small-gain limit. However, it is more interesting for application is to explore 
large-gain regime of ICSRL that needs further investigations. 
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